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ABSTRACT 
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processes increases at larger momentum transfer. 
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1 Introduction 



It seems quite natural to expect that the harder a high-energy collision is, the higher is the 
number of fragments. One of the most tractable and widely explored processes where this 
phenomenon can be seen is deeply inelastic scattering (DIS) with a possibility to change a hard 
scale (Q 2 , gauge boson virtuality) and to detect its influence (if any) on the hadronic invariant 
mass (W) "efficiency". The simplest measure of this efficiency is the multiplicity of secondaries. 

In 1969 Yang and his collaborators basing themselves on the "fragmentation picture" of 
violent collisions, made a qualitative prediction: "...for larger values of the momentum transfer 
t, the breakup process favors larger multiplicities of hadrons" (at fixed hadronic mass). Early 
searches for this effect were inconclusive in both theory and experiment 0. 

In the framework of QCD a quantitative result has been obtained in Refs. @, |j: it appears 
that QCD gluon bremsstrahlung leads to an increase of the hadron multiplicities in DIS, with 
an increase of Q 2 at fixed hadronic mass W, but that this increase is very slow. The distinctive 
feature of this result is that (n) DIS (W, Q 2 ) has a finite limit at Q 2 — > oo and W fixed. 

Later, another result has been claimed in Ref. f|, which predicted an infinite and quite rapid 
growth of (n) DIS (W, Q 2 ) with Q 2 . In the course of the inference of this result it was supposed 
that the influence of the (non-perturbative) composite structure of the nucleon is negligible 
while in []3], |J it plays a key role in the slowness of the Q 2 dependence of the multiplicity at 
fixed W. There is even more trivial objection. On general grounds, the infinite growth with 
Q 2 at fixed W is impossible because of the apriori kinematical bound 

(n) DIS (W,Q 2 )<^ (1) 
m 

where m is some effective mass. 

Quite recently, a weak dependence on Q 2 in the framework of the Dual Parton Model was 
mentioned in ||. 

Experimentally a statistically significant effect of the slow growth of (n) DIS (W, Q 2 ) was 
established in v{y)p interactions and in fi + p interactions ||. The results of the EMC || 
have been described in the framework of QCD in Ref. ||. 

However, subsequent measurements at HERA (HI) were interpreted as a practical Q 2 in- 
dependence [|l(J of (n) DIS (W, Q 2 ) (for the current hemisphere in the hadronic c.m.s.), while 



HI JT21 [T3|] and ZEUS |TJL [T5j reported quite fast Q 2 dependence for the current hemisphere in 
the Breit frame. It should be noted, however, that this last result concerns different bins in W 
for changing Q 2 values. Anyway, the situation is controversial and therefore very interesting. 

In this paper we give our own interpretation of the HERA data on charged hadron mul- 
tiplicities in the current fragmentation region; as will be seen in the text below, these are in 
agreement with Yang's general hypothesis Jl| and our early QCD results f||, f|] (see also the 
review |TT 
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2 Hadronic Spectrum and Multiplicity in DIS 



According to the factorization for inclusive spectra in DIS, the hadronic spectrum in DIS is 
represented by two terms: 

dn . 2 \ f l dz ^2\dn , TJT ^ n o 



dy 



{W ^ ) = l^ w{x , z , Q Y^W^ ) + ^, (2) 



where x = x + (1 — x){rrih/W) exp(— y), y is the rapidity of the detected hadron and m h is 
its mass. In Eq. (Q) dn/dy defines the hadronic spectrum in partonic subprocess, while the 
quantity dn^/dy describes the spectrum of the proton remnant. The latter does not contribute 
to the current fragmentation region of DIS at HERA energies. 

Correspondingly, the average hadronic multiplicity in DIS is represented by 

(n) DIS (W,Q 2 ) = / -w(x,z,Q 2 )(n)(W eff ,Q 2 ) + (no). (3) 

J xq Z 

For small x the weight w(x, z, Q 2 ) in (0), (^) is of the form: 

x 



w(x,z,Q 2 ) = D«[-Q^Q 2 )f g (z,Qi) 



z 



x I J^Dl{^,Q 2 ,Ql)f 9 (z,Q 2 )\ \ (4) 



As can be seen, the hadronic spectrum in partonic subprocess, dn/dy, and the hadronic multi- 
plicity (n) depend on the effective energy, which is smaller than W: 



Kf= Z T^ W2 - ( 5 ) 



In what follows, we shall work in the c.m.s. of the final hadrons. In terms of rapidity, the 
current region in the c.m.s. corresponds to 

- Y < y < (6) 

(it is assumed that the proton goes in the positive direction). 

In our papers P, [| it has been established that the total hadronic multiplicity in the partonic 
subprocess of DIS is related to the hadronic multiplicity in e + e~ annihilation: 

(h)(W,Q 2 )^( n y +e -(W) (7) 

(up to small NLO corrections, which decrease in Q 2 ). 

In the partonic subprocess, the rapidity varies in the range 

-Y<y + y <Y 7 (8) 
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where Y = \n(W e ff/m h ) and 

yo = ~in(^). O) 



2 \l-z. 

The quantity yo determines the rapidity of the centre of mass of the partonic subprocess in the 
centre of mass of the complete process. On integration over z, the region (^) is "smeared" into 
the region 

- Y < y < Y, (10) 

with Y = ln(W/m h ). 

The average value of the effective energy in (|^), (|3]), available for particle production, ap- 
peared to be dependent on both W and Q 2 ^ |J: 

(W eff ) 2 ~ k{Q 2 )W 2 . (11) 

The efficiency factor k(Q 2 ), which stands in front of W 2 in ([H]), is much less than 1 and grows 
slowly in Q 2 . 

From formulas ([3D, (0) and (|TT|) , one can see that the rise of the average hadronic multiplicity 
in DIS has the same physical nature as in e + e~ annihilation. For the first time this behaviour 
has been experimentally established by HI in 1996 |10| . 



However, the QCD growth of (n) is delayed in DIS by the bound-state effects and the slow 
QCD evolution of the structure function. This is why we predicted that the Q 2 dependence of 



(n) at fixed W should remain numerically weak at HERA energies |9|, |TT 



It follows from (|]), (|9|) that the centre of the spectrum is shifted to the region of positive 
rapidities and tends to zero at asymptotically high Q 2 M: 



(yo) 



Q 2 ->oo 



ln(hiQ 2 ))' 



(12) 



The hadronic spectrum in partonic subprocess in the c.m.s. of DIS has the form 

dii^ 1 4- - — 

— = n^-(W eff ,Q 2 )D h (W eff ,y). (13) 
Normalization in the RHS of Eq. (|TJ) is done in agreement with formula (^). 

3 Hadronic Multiplicities in the Current 
Fragmentation Region 

To calculate the multiplicity of charged hadrons in the current fragmentation region, we have to 
define expressions of the quark distribution at small x, of the hadronic spectrum in the partonic 
subprocess as well as of the multiplicity of charged hadrons in e + e~ annihilation. 
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For the quark distribution, we use an analytical expression from Ref. ||16|| , in the case of 
soft initial conditions. As was shown in flEfl , at small x it is in good agreement with the data 
on the structure function from HERA in the wide ranges of Q 2 . Namely, at high Q 2 , we have 



^(z,Q 2 )~r/ 1 (t)exp(-^/2), 



with d = & + 2QN f /27. The variable 



t = 2W6£ln 



is related to the QCD evolution parameter 



(14) 



(15) 



A> WW 

where /3q = 1 1 — 2 Nf /3 is the /3-function in lowest order and 

t 

T ~ 21n(l/z)' 



(16) 



(17) 



The quark and gluon distributions from Ref. |T6| obey the GLAP evolution equations [17 . 



The expression of the initial gluon distribution at z closed to 1 is chosen to have the following 
form 



fo{z, Qq 



<0) |2— 1 

We have omitted constant factors in the RHS of Eqs. fll4|) and (|H|) as they do not influence 
our final results. 

The spectrum of hadrons in the partonic process D h was calculated by many authors. We 
use the expression from Refs. |l8j (N is a normalization factor): 



D\WX) 



calculated in the variable 



N 



exp 



-k - -s5 - -(2 + k)5 2 + -sd 3 + —kd 4 
8 2 4 V ' 6 24 



E h 



(19) 



(20) 



Here Eh is the energy of the detected hadron. 

The average value of (, (q, and its dispersion a are given by the formulas: 



Co 



-T 1 + 



P 48 



24V/3or^ 



1 &rJ' 



a 



3 I 48 



1/4 



1 _ Po I ^ 



64 V 



1 + - 

8r 



(21) 
(22) 
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where 



and 



r = In ( ^ ) (23) 



p [3 ( 48 \ 1/4 / u \ 



16U\PorJ V 1 + 47j' ^ 

(25) 

(26) 

u 

Here p = 11 + 2JV//27, = 1 + JV)/27. 

At low (effective) energies we use the fit of the low-energy data on multiplicity of charged 




hadrons in e + e annihilation from Ref. 19 



(n) ee = 2.67 + 0.48 In W\ (27) 



while for high energies (W e ff > 10 GeV) we apply the fit from Ref. |20|, which well decsribes 
e + e~ data up to LEP energies: 

(n) e+e ~ = -1.66 + 0.866 exp(1.047v / ln W 2 ). (28) 

We have corrected ( p7|) for a fraction of the charged particles from K® and A (A) decays. 

Figure 1 represents the result of calculations of charged multiplicity in current hemisphere 
of the c.m.s. by using formulas (0) and ([L4]) (solid curves) in comparison with the HI data 
from Ref. [ 10 1 . As can be seen, our QCD predictions are in very good agreement with the data. 
These are quite compatible with a slow growth of (n) DI (W, Q 2 ) in Q 2 at fixed W. 

Figure 2 demonstrates a rapid rise of {n) DIS (W, Q 2 ) in the variable W for different values of 
Q 2 , which was predicted many years ago in Refs. and seen previously in e + e~ annihilation 
(the very values of Q 2 taken from [0). Let us note that the HI data presented in Fig. 2 (see 
Table 4 in [|IIJ]) do not correspond to some fixed values of Q 2 , in contrast with the experimental 
points in Fig. 1. 

In order to obtain multiplicity of charged hadrons in current region of the Breit frame, the 
c.m.s. spectrum ([□!) must be integrated in the region 

- Y < y < y B , (29) 

where 



The quantity 



4x(l-x) (31) 
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in the RHS of Eq. (^) is the velocity of the Breit frame in the c.m.s. So, ys corresponds to 
zero rapidity in this frame. 

The results of our calculations of the multiplicity of charged hadrons in the current region of 
the Breit frame are presented in Fig. 3 as a function of Q 2 , in comparison with the HI data [R| 
(solid squares) and ZEUS data |T5j (solid circles). 

The theoretical curves in Figs. 1-3 correspond to the following values of parameters: 

Ql = 1 GeV 2 , A = 0.25 GeV. (32) 
The parameter n g = 6.1 in Eq. (|I9| ) is taken from one of the MRST sets of parton distribu- 



tions |21| 



It should be noted that the strong Q 2 dependence seen by HI and ZEUS in the Breit frame 
has nothing to do with the Q 2 dependence of (n) DIS (W, Q 2 ) in the c.m.s. (see Fig. 1); to a large 
extent it has a kinematical origin. The point is that an increase of Q 2 at fixed W is equivalent 
to an increase of x. As a result, the current region of the Breit frame (|29| ) enlarges. Thus, the 
rapid growth of hadronic multiplicity in the Breit frame in Q 2 (at fixed W) reflects a strong 
increase of the hadronic spectrum towards the central region. 

As for the increase of (n) D (x, Q 2 ) in Q 2 at fixed x in the Breit frame, it has been found that 
it is similar to that in e + e~ annihilation at high Q 2 , while there is a discrepancy between DIS 
and e+e" data at low Q 2 P]-flI5[. Within our approach, it can be understood as follows. On 



the one hand, the increase of Q 2 results in an increase of the height of the spectrum (because W 
grows). On the other hand, the position of the spectrum, (y ) as defined in (g), tends towards 
the region of positive rapidities. 

These two phenomena go in opposite directions. At Q 2 high enough, (yo) varies very slowly 



with Q 2 (12). As a result, the rise of hadronic multiplicity in the Breit frame is analogous to 
that in e + e~ annihilation. At low Q 2 , (yo) changes more significantly |4[]. This effect partially 
compensates the growth of the spectrum and there appears a significant difference between DIS 
and e + e~ data. 

Finally, it is interesting to analyse the case when x increases while Q 2 remains fixed. In this, 
W decreases, which results in a rapid decrease of the spectrum. At the same time, however, 
the current region in the Breit frame becomes larger in accordance with formulas (|2"§|), flUUp. 
The effects are of the same order but opposite in sign. The ZEUS data in the Breit frame (see 
Table 2 in Ref. |D|) show that there is a slow rise of the hadronic multiplicity in the current 
hemisphere in the variable x at different fixed values of Q 2 . 



4 Conclusions 



In this paper we have presented the results of QCD calculations of the multiplicity of charged 
hadrons in the current hemisphere of DIS. Both the c.m.s. and the Breit frame are considered. 
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We have shown that the HI data are in agreement with Yang's hypothesis and our QCD 
predictions. Namely, the efficiency of high energy collisions does weakly depend on the hard 
scale of the process (momentum transfer Q 2 ). It means that at fixed energy the efficiency of 
a particle production in hard processes increases with the shrinking of the interaction region 
(~ 1/Q, in DIS). 

The observed rise of the hadronic multiplicity with Q 2 in the current region of the Breit 
frame has both a dynamical and kinematical origin. This is why a direct comparison of available 
DIS data in the Breit frame with e + e~ data is not completely correct. 

New data from HERA on the hadronic multiplicity in the current region of the c.m.s. as 
well as measurements of the total multiplicity in DIS as functions of two variables (Q 2 and 
W/x) would be very important. 
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Figure 1: The Q 2 dependence of the multiplicity of charged hadrons in the current fragmentation 
region of the c.m.s. in intervals of W 
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Figure 2: The W dependence of the multiplicity of charged hadrons in the current fragmentation 
region of the c.m.s. at fixed values of Q 2 



10 




2 3 4 

10 10 10 10 

Q 2 (GeV 2 ) 



Figure 3: The Q 2 dependence of the multiplicity of charged hadrons in the current fragmentation 
region of the Breit frame at fixed values of W 
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